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ABSTRACT 


An  elastic  wedge  of  interior  angle  x IT  is  subjected  to  spatially 
uniform  but  time  dependent  shear  tractions,  which  are  applied  to  one  or 
both  faces  of  the  wedge,  parallel  to  the  line  of  intersection  of  the  faces. 

The  transient  wave  propagation  problem  is  solved  by  taking  advantage  of 
the  dynamic  similarity  which  characterizes  problems  without  a  fundamental 
length  in  the  geometry.  The  shear  stress  Tqz  is  evaluated,  and  it  is 
found  that  the  singularity  near  the  vertex  of  the  wedge  is  of  the  form 
r ( 1  /x)  The  results  show  that  the  stress  is  not  singular  for 

interior  angles  less  than  n  .  Asa  special  case  we  obtain  the  dynamic 
shear  stress  generated  by  the  sudden  opening  of  a  semi-infinite  crack  in 
a  homogeneously  sheared  unbounded  medium. 


The  results  of  this  paper  were  obtained  in  the  course  of  research 
sponsored  under  Contract  No.  N00014-67-A-0109-0003,  Task  NR 
064-496  by  the  Office  of  Naval  Research,  Washington,  D.  C. 
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INTRODUCTION 


Problems  that  are  concerned  with  the  propagation  of  small  de¬ 
formations  in  linearly  elastic  solids  are  generally  solved  by  means  of 
Fourier  transform  techniques.  If  the  region  in  space  is  unbounded,  and  if 
no  characteristic  length  of  the  geometry  enters  the  formulation,  it  may 
reasonably  well  be  expected  that  a  closed-form  solution  can  be  worked  out. 

In  obtaining  this  solution  the  use  of  Fourier  transforms  becomes,  however, 
less  attractive  for  more  complicated  regions  in  space,  such  as  wedges.  In 
this  paper  we  consider,  therefore,  an  alternative  method  of  solution  which  is 
based  on  the  dynamic  similarity  that  characterizes  problems  without  a  funda- 

O 

mental  length.  This  method,  which  has  been  used  extensively  in  supersonic 
aerodynamics,  [1],  [2]  ,  was  applied  by  Miles  to  wave  propagation  problems 
in  homogeneous  elastic  solids  [3]  . 

We  consider  the  transient  waves  generated  by  spatially  uniform 
but  time  dependent  shear  tractions  which  are  applied  to  one  face  of  a  wedge 
in  a  direction  parallel  to  the  line  of  intersection  of  the  faces.  The  other  face 
is  assumed  free  of  surface  tractions.  Once  the  solution  to  this  problem  has 
been  obtained  we  can,  for  arbitrary  vertex  angles,  easily  construct  the 
solution  for  the  cases  where  both  faces  are  subjected  to  shear  tractions,  or 
where  one  face  is  clamped.  As  an  interesting  special  solution  we  obtain  the 
dynamic  shear  stress  generated  by  the  sudden  opening  of  a  semi- infinite 
shear  crack  in  an  unbounded  medium.  For  the  general  case  special  attention 
is  devoted  to  the  singularity  of  the  shear  stress  near  the  vertex  of  the  wedge. 
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The  transient  diffraction  of  plane  waves  by  a  wedge  in  an  acoustic 
medium  was  treated  earlier  by  essentially  the  same  method  by  Miles  [4] 
and  Keller  and  Blank  [5]  . 
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I 


GOVERNING  EQUATIONS 


A  homogeneous,  isotropic,  elastic  wedge  of  interior  angle  xff, 
see  Fig.  1,  whose  faces  are  defined  by  0  =  0  and  6  =  Xlt  ,  respectively, 
is  subjected  on  the  face  0  =  0  to  a  uniform  but  time  dependent  shear 
traction  Tgz  .  For  the  time  being  we  assume  x  s  ^  ;  the  case  x  <  £  will 
be  considered  later.  The  shear  traction  generates  horizontally  polarized 
shear  motion  in  the  z-direction  only,  which  is  governed  by  the  equation 

IS.  Sw .  1  S^w  _  1  S^w  . 

"r  *Sr  r  "dr  '  +  ~2  1  ~  T  — T  »  ( 

r  S  0  c  St 

where  w  is  the  displacement  in  the  z-direction,  c  =  (/i/p)^  is  the  velocity 
of  shear  waves,  and  r  ,  0  ,  z  are  cylindrical  coordinates.  It  is  assumed 
that  the  wedge  is  at  rest  prior  to  t  =  0 

t  <  0  w{r,  0,  t)  =  w(r,  0,  t)  s  0  .  ( 

The  displacement  field  generated  by  a  uniform  surface  traction  of  arbitrary 
time  dependence  can  be  obtained  by  Duhamel  superposition,  once  the  dis¬ 
placements  for  a  surface  traction  varying  with  time  as  the  Dirac  delta 
function  have  been  found.  We  thus  first  consider  the  following  boundary 


cond  itions 


0=0,  r  £  0: 


0  =  xir ,  ri  0: 


_  £  dw 

T0z  r  T§ 


T  =  * 

T0z  r  1T§ 


=  Tj  6(t) 
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The  problem  at  hand  consists  of  finding  a  solution  of  Eq.  (1),  satisfying 
the  initial  conditions  (2)  and  the  boundary  conditions  (3)  and  (4). 

Some  observations  on  the  pattern  of  waves  propagating  in  the 
wedge  can  be  deduced  from  elementary  principles  of  wave  propagation. 

The  surface  traction  (3)  generates  a  plane  wave  with  constant  displacement 


w 


1 


(5) 


This  wave  is  called  the  primary  wave,  and  in  Fig.  1  its  wavefront  at  an 
arbitrary  time  t  is  indicated  by  BD.  Since  the  wedge  is  at  rest  prior 
to  time  t  =  0  ,  the  medium  is  undisturbed  ahead  of  the  wavefront  BD,  and 
as  discussed  above,  the  displacement  is  constant  behind  it.  In  addition  to 
the  primary  wave,  the  vertex  of  the  wedge,  as  well  as  the  non- uniformity 
of  the  surface  traction  across  the  vertex,  generates  a  cylindrical  wave 
with  center  at  O  .  Since  the  displacement  is  continuous  across  the  cy¬ 
lindrical  wavefronts  we  have  for  x  2  ^  :  w  ^  0  along  BC. 

There  is  no  characteristic  length  in  the  problem,  and  it  is  thus 
to  be  expected  that  the  solution  shows  dynamic  similarity,  i.  e.  ,  the  dis¬ 
placement  is  a  function  of  r/t  and  0  .  It  is  then  expedient  to  introduce 
as  a  new  independent  variable 

s  =  r/t  .  (6) 

Asa  function  of  s  and  0  the  displacement  w(s ,  0)  satisfies 


,  2  .2  ,  2  ,  ,2 
Zn  s  .  d  w  ,.  2 s  .  dw  d  w  _  n 

s  (1  -  r-z  +  5(1  -  — r»  jj +  — z  - 0 

c  d  s  c  &  0 


(7) 
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For  s  <  c  ,  Eq.  (7)  is  elliptic.  Upon  introducing  Chaplygin's  transformation 


8  =  -  cosh’  1  (■§) 

S 


(8) 


Eq.  (7)  reduces  to 


(9) 


The  solution  of  Laplace's  equation  may  be  written  as  the  real  part  of  an 
analytic  function  x(8.0)  » 

w(fl,  0)  =  Re  x(8. 0)  .  (10) 

For  s  >  c  ,  Eq.  (7)  is  hyperbolic,  and  may  be  reduced  to  the  canonical 
form  by  the  transformation  s  =  c  sec  Of  . 

From  Fig.  1  the  region  in  which  Eq.  (7)  is  elliptic  is  now  identified 
as  the  cylindrical  region  ABC.  All  that  now  remains  to  be  done  is  to  find  a 
harmonic  function  w(9 , 0)  in  the  segment  0  s  M  xff  ,  s  £  c  satisfying 
boundary  conditions  which  for  take  the  form 

0=0,  s^c:  -|-g-  -  0  (11) 


In  Eq. 
zero, 


c  T 

s  =  c,Oi0s|:  w  =  Wj  =  -  (12) 

s  =  c,^i0s:xtr:w  =  O  (13) 

0  =  Xff  ,  s  i  c  :  =  0  .  (14) 

(11)  we  have  used  that  for  t  >  0  the  surface  traction  has  returned  to 
see  Eq.  (3).  The  function  w(0,0)  can  be  obtained  in  several  ways;  here 
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we  elect  to  map  the  segment  s  £  c  ,  0  £  0  £  xff  on  the  half-plane  T)£  0 
by  a  conformal  mapping  which  was  introduced  by  Craggs  [2] 

C  =  £  +  iT)=  sech  [(3+i0)/x]  .  (15) 

Equation  (15)  can  also  be  written  as 

£  =  [cosh  —  cos  —  +  i  sinh  —  sin  —  ]  .  (16) 

X  X  X  X 

The  mapping  of  the  segment  s£c,0  £0£xff  on  the  £  plane  is  shown 
in  Fig.  2,  where  the  positions  of  the  various  points  are  indicated.  The 
boundary  conditions  ( 1 1 )- ( 1 4)  are  converted  into  conditions  on  the  real  axis, 
and  we  find  for  x  s  1 


U 1 

<  1  : 

O 

II 

(17) 

1  £ 

i  *  l/cosfx  : 

W  =  Wj 

(18) 

1/cos  : 

2  X 

w  =  0 

(19) 

i  * 

-  1  : 

w  =  0 

(20) 

For  1  &  x  ^  i  the  point  B*  is  located  on  the  negative  real  axis,  and  (17)- 
(20)  must  be  modified  accordingly. 

In  this  paper  we  are  interested  in  the  stresses,  in  particular  in 
the  shear  stress  Tgz  »  which  may  be  written  in  the  form 

T0*  *  f  Re  If  ]  <21> 

According  to  (17)  d  \  /d£  is  real  for  |$  |  s  1  ,  and  from  (18)-(20)  we  con¬ 
clude  that  d^/dC  is  imaginary  elsewhere  along  the  real  axis,  except  at 
the  point 
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(22) 


?  =  «B  =1/cosfi 

At  this  point  Re  \  is  discontinuous,  and  Re  d\/d£  is  a  delta  function, 
which  implies  that  d)(/d£  has  a  simple  pole  at  £  =  .  An  expression 

satisfying  the  foregoing  requirements  is 


B 


(23) 


<C  -1) 


B 


The  constant  B  is  found  by  integrating  counterclockwise  along  a  small 
semi-circle  around  £  =  ,  and  equating  the  results  to  -  Wj  .  We  obtain 

cr, 


w 


1  ,*2 


B  =  —  -  1)*«  -  -  «...  , 

IT  SB  '  TT  M  2  X 


1  tan  ■?- 


(24) 


From  Eq.  (15)  we  compute 


tanh  (£  +  i  £) 
X  X 

cosh  (£■  +  i  — ) 
X  X 


(25) 


Upon  substitution  of  (16)  and  (24)  into  (23),  and  subsequent  substitution  of 
(23)  and  (25)  into  (21)  we  obtain 


0z 


.  CTi  sin<&> 


Se  (i/[  cos  -  coshfj  +  i  |)]) 


IT  X 


(26) 


The  real  part  can  easily  be  evaluated,  and  the  result  is 


where 


CTi  8in<r;>  r 

Tez  =  -—  -rr*-Fi£’e> 


sinh  —  sin  — 


p  I r  Q\  -  x  x  _ 

l'ct*  '  .  i  T1  e A  .  .  ,8  572 

(cos  •=—  -  cosh  —  cos  — )  +  (sinh  —  sin  — ) 
c  X  XX  XX 


(27) 


•  (2  8) 
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Equation  (27)  represents  the  stress  within  the  cylindrical  wavefront,  r  £  ct  . 

For  r  s  ct  and  0  s  0  s  |  we  have 

T  =  T,  cos  9  6(t  -  —  sin  0)  .  (29) 

re  l  c 

It  is  of  particular  interest  to  investigate  for  t  >  0  the  singularity 

of  as  r  -•  0  .  To  this  end  we  evaluate  asymptotic  expressions  for 

Oz 

sinh(0/x)  and  cosh(8/x)  for  small  .  Using  a  well-known  representation 
for  cosh"  *  ,  we  find  from  (8)  and  (6): 


8=  -  in  {^+  [(ff  -  1]*)  .  (30) 

Equation  (30)  is  subsequently  used  to  write 

I  2  .11  I 

cosh  !  =  H(^)*{l  +  [l-(^)  ]*)  *  +  <-^)*  {l  +  [l-(^)2]t}X  ]  (31) 

i  .11  1 

sinl.  f  =  i[(^)X  (l  +  [l-(^)2]4)’X-(^)’1  (l  +  [l-(^)2)ilX  ]  •  (32) 

Thus  for  (4>«  1 
'ct 


cosh  — 


,i  i. , 

.ct.x  2X 
'  r  ' 


sinh  — 


-  -  -  1 
(ct.x  2x 

'  r  ' 


(33) 


(34) 


For  x  a  ^  the  singularity  in  the  shear  stress  is  thus  obtained  as 


T0z 


2  c  T, 


ff  X 


(2ct) 


l 

x 


sin(^) 


“"<!> 


(35) 
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We  note  that  the  singularity  vanishes  for  0=0  and  0  =  x  ff  »  and  reaches 
a  maximum  for  0  =  xtr/2  .  The  singularity  disappears  altogether  if  x  ^  1  ■ 
The  shear  stress  is  thus  singular  only  if  the  interior  angle  of  the  wedge  ex¬ 
ceeds  it  . 
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BOTH  FACES  SUBJECTED  TO  SHEAR  TRACTIONS 


If  the  faces  defined  by  0  =  0  and  0  =  xff  are  subjected  to  shear 
tractions  Tq?  =  T  ^  6(t)  and  Tqz  =  T 2  6(t)  ,  respectively,  the  shear  stress 
T«  is  obtained  by  simple  superposition.  We  obtain  for  r  s  ct 


0z 


sin  (fi> 


It  X 


trl  Fl<H’6>  +  T2F2<£-6» 


(36) 


where  Fj(r/ct,  G)  is  defined  by  (28),  and 


F2<ct’  " 


sinh  ^  sin  ® 

~I  tr  7  r~S  6  ,  2  '  8  '  O 

(cos  +  cosh  —  cos  — )  +  (sinh  —  sin  — ) 

^  X  XX  XX 


(37) 


In  addition,  we  have  for  r  s  ct  and  0  s  8  «  ir/2  the  plane  wave  (29), 
and  for  x*T£  0  2  xtr  -  ^  we  have 


Tgz  =  r2  cos  (xtr  -  0)  6 Ct  -  ^  sin(  xff-0)] 


(38) 


It  is  clear  that  for  this  case  the  singularity  is  of  the  form 


2c(ti  +  t2)  -7  -  s  7-1 

T“- - TH -  (2ct)  sin(I^>  sin(x}  r 


0z 


(39) 


Thus,  the  singularity  vanishes  if  the  two  shear  tractions  are  of  the  same 
magnitude  and  sense,  i.  e.  ,  if  Tj  =  -  T  2  • 
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THE  SHEAR  STRESS  FOR  STEP  SURFACE  TRACTIONS 


It  is  of  considerable  physical  interest  to  compute  the  dynamic 
stresses  for  the  case  that  the  surface  traction  varies  as  a  Heaviside  step 
function, 


e  "  0  '  r  2  0  :  T92  =  T,  1(.) 
e  =x  r  2  0  :  Tgz  =  0 

For  Os  s  ^  the  shear  stress  Tgz  is  now  obtained  as  the  sum  of  the 
integrals  of  (27)  and  (29) 


(40) 

(41) 


T qz  -  Tj  cos  01(t--^  sin  0) 


cT,  Bin(f-) 


IT  X 


r^^.eidv  .  (42) 


£ 

C 


For  0  2  |  we  obtain 

c  T .  s in(  y—  )  t 

—  -nr  X  r,<£.e)dv  .  (43) 

r 
c 

The  behavior  of  Tgz  ^or  small  values  of  r  is  found  by  intro¬ 
ducing  the  new  variable 


P 


c  v 
r 


The  integrals  in  (42)  and  (43)  then  become 


sin(Ti> 

n  x 


rt  Ct/r 

1(T  -  l)  !  Fi(P»  6  )  dp 

1  1 


(44) 


(45) 
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We  are  interested  in  approximating  (45)  for  large  values  of  ct/r  .  For 
large  p  the  cosh[g(p)/x]  and  sinh[B(p)/x]  may  be  approximated  by, 
see  (31)  and  (32), 


cosh  |  ~  pX  2X  (46) 

A  i  I-  1 

sinh  B.  „ .  px  2X  .  (47) 

The  integral  can  now  be  evaluated,  and  for  x  a  »  x  t  1  ,  and  r/ct  <<  1 
we  find 

T8z~iiTT)<2ct,’*tlsin(fxlsin'®»  r*'‘  •  <48> 

For  v  =  I  and  r/ct  <<  1  we  obtain 

Tez  Ti  Bin9‘n(f>  •  <49> 

We  note  that  no  singularity  occurs  if  the  interior  angle  of  the  wedge  is  less 
than  if  . 


SPECIAL  CASES 


The  integrals  in  (42)  and  (43)  can  be  evaluated  with  particular 
ease  for  three  special  values  of  x  :  x  =  0.  5  ,  x  =  1  and  x  =  2  . 

For  x  =  0.  5  we  have  a  quarter- space  subjected  to  uniform 
surface  tractions.  From  Eqs.  (27)  and  (42)  it  is  immediately  noted  that 
the  cylindrical  wave  vanishes,  as  it  should,  and  we  thus  are  just  left 
with  the  plane  primary  wave. 

The  case  x  =  1  is  concerned  with  a  half- space.  We  assume 
uniform  surface  tractions  that  are  of  different  magnitudes  for  8-0  and 
8  =  it  ,  respectively.  Equation  (36)  then  reduces  to 


8z 


1  ^T1  +T2^  sinh  B  sin  ® 

r  cosh^  8  -  sin^  0 


(50) 


Equations  (31)  and  (32)  become 


U  o  Ct 

cosh  8  -  — 
r 


sinh  B  =  -  [(“)2  -  1]^ 


(51) 

(52) 


Thus  for  rs  ct: 


(T1  +  T2)(sin  8)[(ct)2  -  r2]^ 
[(ct)2  -  r^  sin^0] 


(53) 


The  cylindrical  wave  vanishes,  of  course,  altogether  for  Tj  =  -  T2  •  For 
r/ct  <<  1  we  recover  the  behavior  shown  in  Eq.  (39)  for  x  =  1  . 
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For  step  surface  tractions  the  integration  of  (36)  yields 


t  ,  2  , 2  2.i  .2.2  2.4 

hi  -  tb-f/Je > ds 


where  we  have  introduced  the  slowness  b  , 


b  =  - 
c 


This  integral  can  be  evaluated  to  yield  for  br  s  t 


(54) 


(55) 


T  +T 
A1  l2 


0z 


{  sin  6  £ 


i  -  .  -1  rt  sin  0  -  br 

-•cos  0sin  [— - r - : — si 

E  Lt -  br  sin  0  J 


i  A  .  -  1  r  t  sm  0  +  br  , 
-f  cos  0sin  [ .  ■ '  "  ■. — x]j 

Lt  +  br  sin  0  J  3 


(56) 


For  0  £  ^  the  shear  stress  is  then  obtained  as 


Tgz  =  Tj  cos  0  l(t-br  sin0)  +  I0z  Kt-br) 


(57) 


For  0  i  ~  we  have 

Tgz  =  -  T2  cos  0  l(t-br  sin0)  +  lgz  l(t-br)  .  (58) 

For  br/t  <<1  the  behavior  is  as  indicated  by  (49). 

The  case  x  =  2  corresponds  to  a  semi- infinite  slit  in  an  un¬ 
bounded  medium  subjected  to  different  shear  tractions  on  the  two  faces.  For 

=  T2  =  T  ,  (59) 


Eq.  (36)  reduces  to 


14 


2^ 


T  sinh  9  sin  6  (1  +  cosh  8  +  cos  0) 
(cosh  8 )  -  (sin  0  ) 


For  x  =  2  ,  Eq.  (32)  becomes 


(60) 


sinh  $  0  =  -  2"  ^  (-y  -  1)^ 


(61) 


After  some  further  manipulation  we  then  obtain  for  r  <  ct 


T,-t.ct  ,C°St(f-6)  Sin^(f~  W 

T0  z  tr  r  *■  t-(r/c)  sin  0  "  t+(r/c)  sin0 


} 


(62) 


The  behavior  for  r/ct  <<1  agrees  with  what  is  obtained  from  Eq.  (39) 
for  x  =  2  . 

The  integrals  to  determine  the  stress  for  step  surface  fractions 
can  again  be  evaluated  explicitly,  and  we  obtain  for  br  £  t 


(T^  IT  =  {(t-br)^t[br(l-sin0)]^tan~1  [brj~^:ne)  ]*}cos$(g- 


6) 


{(t- br)^+  [br(l  +  sin  0)]^  tan"  1  l^YTsin  0)  ^sin£  (63) 


where  T  is  the  applied  surface  traction.  For  0  i  0  i  |  and  2  v  <:  0 

the  shear  stress  Ta  is  obtained  as 

oz 


For  |  s  0  S  we  find 


Tqz  =  T  cos  P  l(t-br  sin0)  +  l(t-br) 


0z 


(64) 


TSz  =  J8z  1<t-br| 


(66) 


If  r/ct  >  v  1  the  singular  behavior  is  of  the  form 


n  t 


2T 


r  ~  (f7>*  {«,*(!- ei  -  Binttf-en 


(66) 
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which  agrees  with  (48). 


Equations  (64)  and  (65)  when  superposed  on  a  uniform  shear  stress 
T  yield  the  solution  for  the  stress  Tfi  due  to  the  sudden  opening 


of  a  semi- infinite  crack  in  an  unbounded  medium  which  was  in  a  state  of 


uniform  shear  prior  to  fracture.  As  is  common  in  crack  problems  we  find 
that  the  stress  is  singular  near  the  crack  tip  as  r"  ^  . 
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CONCLUDING  REMARKS 


Although  only  the  case  x  £  0.  5  was  treated  here,  it  is  clear 
that  wave  propagation  in  wedges  with  interior  angles  less  than  *r/2  can 
be  treated  by  using  superposition  in  conjunction  with  symmetry  or  anti- 

jj{  5J5 

symmetry  properties.  Thus  if  the  face  0  =  x  if  (0.  25  £  x  <  0.  5)  is  free 

of  traction,  and  9=0  is  subjected  to  T ^  6(t)  ,  we  can  simply  use  (36) 

$  $ 

with  x  =  2x  ,  and  to  obtain  the  solution.  If  the  face  0  =  x  IT 

is  clamped,  i.  e.  ,  w  =  0  ,  we  substitute  x  =  2x  and  T-,  =  in  Eq.  (36). 

♦♦  * 

For  x  <  0.  25  the  solution  for  0.  25  £  x  <  0.  5  has  to  be  used  in  the 

just  described  procedure. 

It  is  finally  concluded  that  the  method  of  homogeneous  functions 
in  an  efficient  method  to  study  the  dynamic  response  of  an  elastic  wedge 
to  spatially  uniform  shear  tractions.  For  surface  tractions  varying  in 
time  as  Heaviside  step  functions  it  was  found  that  for  interior  angles  xff  , 
the  singularity  of  the  shear  stress  Tqz  is  of  the  form  [ r^1  ^  "  1 /( 1- x)l, 
which  shows  that  the  stress  is  singular  only  for  x  2  1  . 
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